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Abstract 

The phase structure of lattice QCD with Wilson fermions is discussed. An- 
alytic and numerical evidences are given for the spontaneous breaking of parity 
and flavor symmetry, which naturally explains the existence of the massless 
pion at the critical hopping parameter K c without recourse to the chiral sym- 
metry absent in the Wilson fermion formulation. New numerical evidences 
are also presented for the multiple structure of the critical lines in the weak 
coupling region. A connection between the phase structure and the finite 
temperature phase transition is briefly mentioned. 



*talk presented at "QCD on massively parallel computers" held at Yamagata, Japan, on March 
16-18, 1995. 
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1 Introduction 



Elucidating the nature of chiral transition separating the high-temperature quark- 
gluon plasma phase from the low-temperature hadron phase has been one of the focal 
points of effort in recent numerical simulations of full lattice QCD including light 
dynamical quarks. Much of work in this direction has employed the Kogut-Susskind 
quark action since it retains a U(l) subgroup of chiral symmetry. On the other hand, 
studies with the Wilson quark action are less well developed in spite of the effort 
over the years |], ||, ||. 

A main problem for studies with Wilson fermions originates from the presence of 
explicit chiral symmetry breaking in the action at non-zero lattice spacing, which is 
introduced to avoid the well-known species doubling problem. Clearly, before serious 
studies on the phase structure at finite temperature are initiated, more understanding 
on the effect of the explicit chiral symmetry breaking at zero temperature is necessary. 

Let me consider a phase space of (3 = 6/g 2 and K = 1/2M = l/(2m g a + 8) 
at zero temperature, where g is the bare gauge coupling, m q the bare quark mass, 
and a the lattice spacing. There is much evidence from analytical considerations 
and numerical simulations of hadron masses that the pion mass vanishes along a 
line K = K c (/3), the critical line, which runs from K c (/3 = 0) = 1/4 in the strong- 
coupling limit to K c (/3 = oo) = 1/8 in the weak-coupling limit. However the massless 
pion appeared on K c (f3) can not be regarded as the Nambu-Goldstone boson of the 
spontaneous breaking of chiral symmetry since it is explicitly broken in the Wilson 
fermion formulation at non-zero lattice spacing. 

Why does the massless pion appear then ? It has conventionally been explained 
as follows: The appearance of the massless pion on K c may be an evidence that the 
tuning K effectively recovers the chiral symmetry, explicitly broken by the Wilson 
fermion formulation. Off course, this in NOT an explanation, since it does not tell 
why the massless pion appear on K c . It also can not tell what happens beyond K c . 

Instead of the conventional explanation, a possible explanation was proposed 
some years ago, to explain the appearance of the massless pion[Q, on the basis of 
new phase structure of lattice QCD with Wilson fermions. Later some analytic^]. [| 
and numerical@, |], |§ evidences were given for the existence of the phase structure, 
but it is far from being completely established. In this report, after reviewing the 
previous accomplishment, more evidences are supplied to established the proposed 
phase structure. 

2 Possible explanation 

In this section we review the possible explanation for the massless pion appeared 
on K c in the Wilson fermion formulation^, ||. Let me consider lattice QCD with 
a single flavor first. In the theory there exists one 'V meson", corresponding to a 
pseudo-scalar operator n(x,t) = 'ipij5'ip(x,t). At K < K c the correlation function of 
the operator for large t behaves such that 

£<7r(0,0)7r(x,t)>~Zexp(-t/O (1) 
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where the correlation length £ is related to the mass of the n meson as £ = l/(m T a). 
Since the correlation length diverges on K c the correlation function for large t behaves 
at K = K c such that 

^(7r(0,0)7r(x,t))^Zr a (2) 

X 

with some exponent a. If we go beyond K c (K > K c ), it is expected for large t that 

J2(A0,0)n(x,t)) ~ V(n) 2 + Zex V (-t/0 (3) 

with (7r) 7^ where is a space volume. 

The above situation is summarized as follows. Since the massless tt meson implies 
the divergent correlation length on K c , the operator corresponding to the tt meson 
has the long-range order at K > K c in such a way that (tt) ^ 0. This condensation 
spontaneously breaks the discrete parity at K > K c and the second order phase 
transition from the parity symmetric phase to the parity broken one occurs at K = 
K c . Therefore the massless tt meson can be interpretated as the divergent correlation 
length associated with this second order phase transition, not a Nambu-Goldstone 
boson. 

According to the above interpretation we expect near K c , (m^a) 2 ~ (K c — K) 2u oc 
(M — M c ) 2l/ with another critical exponent v. If we regards the effective theory 
for the tt meson as a four- dimensional scalar field theory, we can expect the phase 
transition to be mean-field like up to the logarithmic corrections and therefore v = 
1/2, which reproduces well-known PCAC relation that {m^a) 2 cx m q a where a quark 
mass m q a ~ M — M c . 

Now we consider QCD with 2 flavors, where members of the pseudo-scalar mesons 
are ir a = tfjij^T""^ (non-singlet) and rj = tpij^ip (singlet). Since the observed mass 
spectra of these mesons satisfy m v > m n , it is expected at K = K c that m n a = 
0,m v a 7^ . To satisfy this, there exists a phase transition at K — K c , which 
implies the condensation such that (7r 3 ) ^ 0, (rj) = at K > K c . The above 
condensation spontaneously breaks both parity and flavor symmetry. In this case 
the neutral 7r° meson corresponding to the 7r 3 operator is identified with the massless 
mode of this phase transition and becomes massless only at the transition point 
K = K c . The charged 7r =1= mesons corresponding to the linear combination of it 1 
and 7r 2 operators are Nambu-Goldstone bosons associated with the flavor symmetry 
breaking: m T o = m n ± in the phase with (7r 3 ) = and m % ± = in the phase with 
(7r 3 ) 7^ 0. On the other hand the 77 meson stays massive for all K since (77) is always 
zero. This phase transition can naturally solve not only the problem of the massless 
pion but also the lattice U(l) problem||. 



One can rigorously prove that (rj) = using a la Vafa-Witten theorem[10|, while 
no such restriction can be applied to the operator 7r 3 , since the Wilson fermion 
determinant is shown to be real-positive for Nf = 2 even in the presence of the 
source term iHijj^/5T 3 ijj. 



3 Evidences for the phase transition 
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3.1 Analytic results 

In the strong coupling (/3 = 0) and large N c limit we can analytically calculated the 
behavior of order parameter and meson masses 0, |5|: 

(-0*75-0) =0 M 2 > M 2 C 

<0z 7 5^> =± 16 _ 2 M2 V / 3(4-M^)x4JV c M 2 <M C 2 

where the critical point is given by M c = ±2(K C = ±1/4) with M = 1/2K. In this 
limit there is no essential difference between Nf = 1 and Nf = 2. It was also shown 
that the ir meson becomes massless on M c ||. ||. 

In the strong coupling and 1/N C expansions the behavior of the order parameter 
is essentially the same for Nf = 1 case: (^r^) ^ at M 2 < M 2 (/3), and the pion 
becomes massless on M c (/3), which was calculated in the expansion of /?|| ||. For 
N f = 2 we obtain (tpi^r 3 ^) ^ at M 2 < M 2 (/5) while {ipij 5 ip} = for all M, 
as expected^, §]. Meson spectrum also agrees with our expectation in the previous 
section. 



3.2 Results from numerical simulations 

To calculate quark propagators on a given configuration iterative methods such as 
Conjugate Gradiant (CG) method are usually used to invert the lattice fermion 
matrix. For the case of the Wilson fermion, if K approaches K c , the inversion 
becomes more and more difficult: a number of CG iteration rapidly increases. To 
overcome this difficulty and to calculate non-zero order parameter, we add an external 
field coupled to the order parameter and take H — > limit after the average over 
configurations is taken. 

Let me consider the quenched QCD first. If we add the external field iH%() r )^il) to 
the Wilson fermion action, it is easy to show that 

Re<#y&V)(tf, H) = -Re(^i l5 ^)(U, -H) (4) 

where (0)(U,H) means a value of an operator O on a given gauge configuration 
U and the external source H. If H — > limit was non-singular, this would lead 
to Re(0i75^) (U, 0) = before the average over gauge configurations. Therefore 
{'ijji^^ip) (H = 0) = in the numerical simulation on a finite lattice, where the limit 
is always non-singular. However this does not means (ipi^^ip) = 0. Since, in the 
parity broken phase, H — > limit is singular as the volume V goes to infinity one 
should take limits in the following order: 

limflim (i)i^il)){H)]. (5) 

Since actual numerical simulations are always performed on finite volumes, one 
should extrapolate to H = from non-zero H, not put if = in the simulation, in 
order to obtain the correct value of (^75 ip) at H = in the infinite volume limit. 
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This external source method has been applied to show that the parity breaking 
phase transition exists at T = 0, where T is temperature, using the quenched QCD 
simulation at /3 — 5.0 on a 4 4 lattice^]. Such a phase transition disappear at T > T c , 
where T c is the critical temperature of the pure gauge theory, using the quenched 
QCD simulation at /3 = 5.6, 5.8 on a 8 3 x 4 lattice 0. 

For the case of full QCD with 2 flavors, a term iHipi^^ip should be added. 
Even in the presence of the term the Wilson fermion determinant is still real-positive. 
Therefore the hybrid Monte Carlo method is applicable to this case. Again we have to 
keep H non-zero for simulations on a finite volume, in order to find the spontaneous 
breakdown of parity and flavor symmetry. 

The hybrid Monte Carlo simulation has been performed for the full QCD with 
N f = 2 at p = 5.0 on 4 4 , 4 2 x 6 2 , and 6 4 lattices^]. Although some evidence for 
the parity and flavor symmetry breaking phase transition was observed, H — > 
extrapolation was too difficult to prove the existence of the phase transition without 
doubt. Lattice sizes used may be too small for this /?. 

Since the numerical evidence of the phase transition is established for the quenched 
case but is ambiguous for the case of the full QCD, we carry out new full QCD simu- 
lation with Nf = 2 at (3 = 0, where the analytic expressions for the order parameters 
are available for non-zero external field H. In Fig. |l|, (tpij^T 3 ^) is plotted as a func- 
tion of M at = on a 4 4 lattice with H = 0.05,0.1,0.2,0.3, together with the 
analytic expressions with the corresponding H in the strong coupling and large N c 
limits, Since the agreement between the numerical data and the analytic predictions 
are very good, it can be established at least in the strong coupling limit that the 
parity and flavor symmetry breaking phase transition really occurs at M — M c = 2 
(K = K c = 1/4) and (^ij^t 3 ^) becomes non-zero in the H — > extrapolation at 
M < M c (K > K c ). We have checked that the finite size effect to the order parame- 
ter is very small by comparing the data on a 4 4 lattice with the data on a 6 4 lattice 
or with the analytic expressions corresponding to the values in the infinite volume 
limit. 

As a further check we have calculated susceptibilities X55, Xoo, an d X50, which are 
defined by 

X55 = ^([E^75T 3 ^(X)] 2 ) (6) 
' x 

Xoo = ^(LE^)] 2 ) (7) 

' x 

X50 = ^([E^7 5 r 3 ^(x)] x EWy)]>, (8) 

* x y 

both analytically and numerically. For example, X55 is plotted in Fig. |2|. The agree- 
ment between them are reasonably good, except some deviations are observed deep 
in the broken phase. Therefore, the result for the susceptibilities also support our 
conclusion about the existence of the phase transition. 

If we consider the above result from the opposite side it can be said that the 
success of our analytic calculations, except the deviation in the susceptibilities deep 
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in the broken phase, establishes the validity of the strong coupling and large N c 
expansions near the limits. Furthermore, since the expansions predict the existence 
of the phase transition, we can conclude that the parity and flavor symmetry breaking 
is shown to exist in the strong coupling region. 

Finally let me comment on the deviations in the susceptibilities observed deep in 
the broken phase. We think that the deviations are caused by the massless Nambu- 
Goldstone charged pions appeared in the broken phase, whose effect to the suscep- 
tibilities can not be included in the analytic calculation of the large N c limit. We 
expect the agreement will become better if we include the effect in our analytic 
expressions. 



4 Phase structure in the weak coupling region 

In the previous section we have established the phase structure in the strong coupling 
region. There are two critical lines M c (/3) and —M c (f3) in f3 - M plane, which separate 
the broken phase from the symmetric phases: (i/jij^t 3 ^) 7^ at — M c (/3) < M < 
M c (/3) and (ipi~f 5 T 3 ip) = at M < —M c {0) or M c ((3) < M. In this section we 
investigate how this phase structure will be modified in the weak coupling region. 

In Fig. [3] we propose the expected phase structure of lattice QCD with 2 flavor 
Wilson fermions in the (3 - M plane. At (3 < (3 W there are only two critical lines, 
while there are ten critical lines at /3 > (3 W and 5 points where two lines meat 
at f3 = oo. Each point corresponds to one continuum limit, whose low energy 
spectra are composed of a part of sixteen zero modes of the lattice fermion: one 
mode near p = (0, 0, 0, 0) appears at M = 4, four modes near p = (tt, 0, 0, 0) and 
its permutation at M = 2, six modes near p = (tt, tt, 0, 0) and its permutation at 
M = 0, four modes near p = (tt, tt, tt, 0) its permutation at M = —2, and one mode 
near p = (tt, tt, tt, tt) at M = —4. Note that the phase diagram is shown to be 
symmetric under M —>■ —M transformation. It is also noted that the width of the 
broken phase in the weak coupling should behave like C(g 2 ) exp[— B/g 2 ] where B 
is some constant and C(g 2 ) is some regular function of g 2 , since the appearance of 
the width has not been detected by the usual perturbative expansion, and therefore 
should be non-perturbative. 

There are several evidences for the existence of the multiple critical lines in the 
weak coupling region. First of all, the phase structure similar to Fig. |3] has been found 
for the Gross-Nevue model in 2 dimensions in the large N limit, regulated on the 
lattice using Wilson fermions 0, |6|, |TH|. The existence of the multiple structure of the 
critical line has been also suggested from studies of eigen- value spectra of the Wilson 



fermion matrix on quenched configurations [fP2| and full QCD configurations [H 

To establish the multiple structure of the critical lines we have performed new 
simulation using quenched configurations without adding the external field iHi/jj^ip. 
We first monitored the quantity so-called "pion norm" : ^l^^)} 2 ) ■ This 

quantity becomes very large in the broken phase, due to the contribution from zero 
modes of the Wilson fermion matrix. Although some signals for the multiple struc- 
ture of the broken phases are observed in this quantity, more clear indicator for the 
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multiple structure of the broken phases is found to be the number of the CG iter- 
ation necessary to invert the Wilson fermion matrix. This quantity becomes very 
large and is proportional to the volume in the presence of zero modes of the Wilson 
fermion matrix, which are expected to exist in the broken phase. We have plotted 
the number of CG iterations as a function of M at (5 — 5.0, 5.5, 5.6, 5.7 on 6 4 lattices 
in Fig.||, and at (3 = 5.8, 5.9, 6.0 on 12 4 lattices in Fig.||. The multiple structure of 
the broken phase is clearly observed at (3 > 5.6. 

Finally we calculate again the order parameter {ipi^^ip) with small non-zero ex- 
ternal field, H = 0.01, in the weak coupling region (/3 = 5.9). As seen from Fig. || 
the multiple structure of the broken phases is again observed 



5 Conclusion and further investigations 

The phase transition which causes or is caused by the massless pion appeared in 
lattice QCD with the Wilson fermion formulation seems to exist at K = K c . In 
the strong coupling region at K > K c it is established that (ipij^ip) ^ for the 
quenched QCD and (ipi^T^ip) ^ for the 2 flavor full QCD. The evidence of the 
multiple structure of the critical lines are also given in the weak coupling region in 
the case of the quenched QCD. 

As mentioned in the introduction, it is very important to study the phase struc- 
ture of lattice QCD at finite temperature. In the case of the Wilson fermion formula- 
tion one should understand, based on the finding of this report, the relation between 
the critical line discussed in this report and the finite temperature phase transition 



line. Recently such an investigation has been carried out[14fl for the 2 flavor QCD. 
The study concludes that the critical line of the massless pion at finite temperature 
forms a cusp at finite gauge coupling, and that the line of thermal transition runs 
past the tip of the cusp without touching the critical line. The detail can be found 
in ref . H . 
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Figure 1: (vpi'y^T 3 ^) as a function of M at (3 = on a 4 4 lattice with if 
0.05,0.1,0.2,0.3. The lines are analytical prediction obtained at (3 = in the lar 
N c limit with the corresponding H. 
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Figure 2: Susceptibility %55 as a function of M at f3 = on a 4 4 lattice with 
if = 0.05,0.1,0.2,0.3. The lines are analytical prediction obtained at (3 = in the 
large iV c limit with the corresponding H . 
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Figure 3: The expected phase structure of lattice QCD with Wilson fermion in (3 - 
M plane. 
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# of CG iterations vs. M 
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Figure 4: The number of CG iterations as a function of M at (3 — 5.0, 5.5, 5.6 and 
5.7 on 6 4 lattices. 
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# of CG iterations vs. M 
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Figure 5: The number of CG iterations as a function of M at /3 = 5.8, 5.9 and 6.0 
on 12 4 lattices. 
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Figure 6: (ipi^ip) as a function of M at (5 — 5.9 on a 12 4 lattice with if = 0.01. 
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